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20  Abstract  (Continued) 


' same  region.  Major  warmings  can  only  develop  when  the  prewarming  lower  stratospheric  winds 
are  strong.  Damping  controls  the  maximum  amplitude  that  a warming  can  attain  and  the  time 
constant  for  its  growth  rate. 

The  evolution  of  m = I and  rn  = 2 warmings  are  very  different.  A m = 1 warming  is 
characterized  by  pronounced  oscillation  of  wave  amplitude  and  mean  flow  that  result  from 
resonantly  trapped  westward  propagating  planetary  waves  moving  in  and  out  of  phase  with  the 
tropospherically  forced  stationary  planetary  wave.  This  oscillation  can  reach  sufficient  amplitude 
to  decelerate  the  zonal  flow  during  a cycle  to  easterlies  and  create  a critical  level.  Although 
formation  of  a mesospheric  critical  level  is  not  required  to  initiate  a warming,  the  development 
and  propagation  of  critical  levels  in  the  middle  stratosphere  atmosphere  is  central  to  the 
evolution  of  sudden  warmings.  During  a m = I event  the  critical  level  forms  initially  in  the 
polar  region  and  advances  equatorward  in  its  development.  But  a m = 2 critical  level  first 
develops  in  the  equatorial  region  and  advances  poleward.  A m « 2 warming  is  also 
cluracterized  by  a sudden  intensification  after  an  initially  slow  growth  in  contrast  to  slowly 
developing  r,-:  = 1 warmings.  Both  m = I and  m » 2 warmings  are  accompanied  by  mesospheric 
cooling  as  a result  of  the  induced  secondary  circulation  in  response  to  eddy  that  transport  to  the 
polar  stratosphere. 

^ Long  term  integrations  with  a steady  forced  m » 1 wave  show  that  the  mean  flow  evolves  to 
a steady,  asymptotic  state  with  net  cooling  in  the  polar  mesosphere  from  planetary  waves. 

But  steady  m = 2 forcing  leads  to  multiple  generation  of  warmings  which  are  similar  to 
stratospheric  vacillations. 
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NUMERICAL  SIMULATION  OF  SUDDEN 
STRATOSPHERIC  WARMINGS 


I.  INTRODUCTION 

The  most  remarkable  dynamical  event  which  affects  the  stratosphere  is  the  sudden  stratos- 
pheric warming  (SSW),  Major  stratospheric  warmings  occur  on  the  average  every  other  year  in 
the  Northern  Hemisphere  (Schoeberl,  1978),  and  during  the  off  years  minor  warmings  are  usu- 
ally observed.  Satellite  observations  indicate  that  warming  events  also  occur  in  the  Southern 
Hemisphere  (Barnett,  1975;  Hartmann,  1977),  but  these  events  do  not  qualify  as  major  warm- 
ings by  WMO  standards.  The  influence  of  the  SSW  is  not  confined  to  the  middle  atmosphere. 
Quiroz  (1977)  and  McGuirk  (1978)  showed  that  the  1976/77  SSW  definitely  affected  tropos- 
pheric weather  patterns.  In  addition  Ramanathan  (1977)  argued  that  the  SSW  may  have  a 
climatic  impact  on  the  radiation  budget  o(  the  polar  regions.  In  view  of  their  meteorological 
significance  a comprehensive  understanding  of  their  origin  and  development  is  clearly  war- 
ranted. 

Planetary  scale  eddies,  the  ultimate  energy  source  for  the  SSW,  are  formed  in  the  tropo- 
sphere through  baroclinic,  orographic,  and  diabatic  processes.  The  most  important  of  these 
eddies  at  mid-latitudes  are  large  scale  quasi-stationary  waves  with  zonal  harmonics  m - 1 and  2 
(van  Loon  et  al.  1973).  During  winter  these  waves  propagate  vertically,  perturb  the  stratos- 
pheric circulation,  and  transport  heat  and  momentum  from  mid-latitudes  into  the  polar  regions. 
Occasionally,  a very  large  planetary  wave  amplitude  pulse  is  observed  in  the  stratosphere,  often 
associated  with  the  development  of  blocking  patterns  in  the  troposphere.  Within  7 to  10  days 
MiinuMrtpi  suhmillcil  March  .10, 
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following  ihe  sudden  increase  in  wave  amplitude  the  temperature  rises  dramatically  over  the 
pole  as  a result  of  the  northward  eddy  heat  flux,  and  the  zonal  winds  reverse  from  westerly  to 
easterly.  Over  single  stations  at  polar  latitudes,  rocket  observations  show  that  the  warming 
appears  fust  at  high  levels  and  then  descends  into  the  lower  stratosphere.  During  very  strong 
sudden  warmings  the  temperature  increase  over  the  polar  region  may  even  penetrate  deep  into 
the  troposphere  (Quiroz,  1977). 

Previous  theoretical  investigations  of  SSW  events  were  reviewed  by  Schoeberl  (1978)  and 
are  briefly  summarized  here.  Matsuno  (1971)  postulated,  on  the  basis  of  his  fairly  realistic 
simulation  of  SSW  events,  that  a sudden  increase  in  planetary  wave  amplitude  in  the  tropo- 
sphere during  winter  triggers  a vertically  propagating  transient  wave  which  decelerates  the  mean 
zonal  flow.  As  the  transient  wave  travels  to  higher  altitudes  and  lower  atmospheric  densities 
the  deceleration  becomes  more  pronounced  and  Anally  easterlies  are  generated.  The  easterly 
wind  region  provides  a critical  level  for  the  stationary  planetary  waves  below  which  the  mean 
zonal  flow  is  further  decelerated.  The  critical  level  descends  in  response  to  the  continuous 
input  of  planetary  wave  energy  from  the  troposphere  until  it  reaches  the  lower  stratosphere 
where  the  SSW  event  becomes  fully  developed. 


In  Matsuno's  explanation  frictional  and  radiative  damping  play  no  importance  role  in  the 
development  of  a SSW  event.  However,  the  absence  of  damping  prevents  the  decelerated 
mean  zonal  flow  from  relaxing  back  to  the  unperturbed  state  after  eddy  forcing  ceases  and 
causes  planetary  waves  to  be  numerically  reflected  from  stationary  critical  levels. 

Geisler  (1974),  using  a one-dimensional  mechanistic  model  of  a SSW  with  realistic  damp- 
ing and  lower  boundary  conditions,  noted  that  the  descent  of  the  critical  level  is  often  ham- 
pered by  distortion  of  the  flow  ahead  of  the  critical  level.  This  distortion  partially  reflects  the 
wave  before  it  reaches  the  critical  level,  reduces  the  energy  input  into  the  critical  region,  and 
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(hus  slows  ihc  descent  of  the  easterly  wind  region  Me  also  found  that  the  wave  amplitude  forc- 
ing at  the  lower  boundary  must  be  sustained  for  at  least  13  days  to  initiate  a warming  in  agree- 
ment with  some  observational  data. 


Molton  (I‘)76)  studied  SSW  events  with  a mechanistic,  primitive  equation  model  that 
included  realistic  values  of  damping  m contrast  to  Matsuno's  {l*)7l)  model  and  obtained  simu- 
lated warmings  which  showed  considerably  different  development.  For  example,  Molton  did 
not  obtain  a double-layer,  m — I warming  and  found  no  evidence  for  a descending  critical 
level  Instead,  after  an  initial  weakening  of  the  mean  zonal  winds  by  upward  propagancy 
switch-on  pulse,  a region  of  easterlies  appeared  simultaneously  throughout  the  30-80  km  polar 
region  Molton  also  noted  a large  deceleration  of  tropical  mean  zonal  winds  appeared  to  be 
associated  with  the  development  of  the  SSW 

From  the  above  summary  and  the  review  by  Schoeberl  (l')78)  a number  important  ques- 
tions emerge  above  the  mechanisms  and  development  of  SSW  events.  For  example,  how 
important  is  critical  level  interaction’’  What  is  the  difference  between  single  and  double  level 
warmings’’  What  is  the  relative  importance  of  forced  stationary  planetary  waves  and  transient 
planetary  waves’.’  Mow  important  is  damping’’  What  is  the  amplitude  of  the  thermal  planetary 
wave  imbedded  in  a warming’’ 


s: 


In  this  paper  we  attempt  to  explore  and  answer  some  of  these  important  questions.  A 
mechanistic  mcHlel  of  the  middle  atmosphere  previously  discussed  in  Schoeberl  and  Strobel 
(1 978a, b)  is  adopted  and  modified  specifically  to  investigate  the  structure  and  development  of 
SSW  events  This  model  is  similar  to  Matsuno's  (1971)  model,  but  with  some  important 
differences  namely,  realistic  frictional  and  radiative  damping  and  a self  consistent  calculation  of 
the  mean  zonal  wind  structure  We  will  show  that  the  most  critical  input  parameters  that 
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mnuence  the  evolution  of  a SSW  event  are  the  magnitude  of  damping  and  the  strength  of  stra- 

I tosphcric  winds  in  the  15-35  km  region.  Formation  of  a me.vosp/uT/i  critical  level  is  not  essential 

1 

for  the  occurrence  of  any  SSW  However  critical  level  formation  is  fundamental  to  the  develop- 
ment of  in  - 2 warmings 
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II.  MODEL 


The  /onally  averaged  cireulalion  of  ihe  middle  aimosphere  (Ih-l.^O  km)  is  eompuicd  wiih 
a guasi-geoslrophie.  numerieal  model  lhai  explicitly  includes  a self  consisicni  calculation  of 
solar  radiative  healing  due  to  ()i  and  ()|  absorption.  Newtonian  cooling.  Rayleigh  friction,  and 
tropopause  boundary  conditions  based  on  climatological  averages  (Schoebcrl  and  Strobel, 
l')’8a.  hereafter  referred  to  as  Paper  I) 

This  circulation  is  governed  by 
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where  the  variables  have  their  customary  meaning.  .Ml  notation  is  defined  in  .\ppendix  A 

The  relevant  equation  for  the  planetary  scale  eddies  and  waves  can  be  derived  by  scaling 
Laplaces’  tidal  equation  (Schoeberl  and  Geller.  1977).  Each  zonal  harmonic  (w)  is  governed  by 
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This  eguaiuin  is  a slinhi  gencrali/aiion  of  ihc  eguaiion  used  b>  Maisuno  (1*171 1 

f or  ihe  /onallv  averaged  How,  [-.iis  (I)  — (5)  may  be  combined  lo  lorm  a single  eguaiion 
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Ciiven  //.  M.  and  boundary  values  for  i and  </> , Eqs.  (5)  and  (b)  form  a complete  sei 
which  may  he  solved  numerically  by  the  method  described  in  Appendix  B The  appropriate 
values  of  and  o lor  the  generation  of  the  mean  zonal  wind  field  are  given  in  Paper  1 Since 
the  emphasis  is  on  the  winter  hemisphere  at  cxtratropical  latitudes,  photochemical  acceleration 
of  the  thermal  relaxation  rate  may  be  neglected  The  quantity  A/ is  included  in  Eqs.  (1)  and 
(b)  to  correct  the  lower  stratospheric  wind  fields  which  are  only  weakly  forced  by  ozone  absorp- 
tion of  insolation  In  the  studies  by  Matsuno  (1971)  and  Holton  (197b)  a two  dimensional 
mean  zonal  wind  model  of  Ihe  middle  atmosphere  was  constructed  from  observational  data  for 
the  initial  conditions  Here  we  calculate  //.  infer  A/  from  observational  data,  and  solve  the 
steady-stale  form  of  Eq.  (b)  without  Ihe  eddy  flux  divergence  terms  to  generate  Ihe  initial  wind 
field 
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f'igure  1 presents  the  principal  zonal  wind  fields  used  in  this  study.  In  Fig.  la  the  wind 
field  is  similar  to  Maisuno’s  (1971)  initial  field  and  is  generated  with  A/ derived  from  observa- 
tional data  and  k,  increased  from  the  values  adopted  in  Paper  I.  When  the  k,  values  from 
Paper  I are  used,  the  wind  field  In  Fig  lb  is  obtained.  If  M is  set  equal  to  zero,  then  the  latter 
wind  field  is  transformed  to  the  result  in  Fig.  Ic.  The  essential  difference  between  the  Matsuno 
(1971)  and  Holton  (1976)  wind  field  (Fig.  la)  and  the  wind  field  in  Fig.  lb  is  the  polar  night  jet 
strength.  The  principal  characteristic  of  Fig.  Ic  is  the  weak  winds  in  the  lower  stratosphere. 

The  appropriate  boundary  conditions  for  Eqs.  (5)  and  (6)  are: 


d(t> 

dr 


6(«)  at  r - 13.2  (108km) 


M = 0 at  « - ± -2^ 

2 

<6  -£■(«)  at  r = 0 (9km) 
rf.  » 0 at  r - 13.2  (108km) 


(7) 


<6  “ 0 at  W “ ± — 

<b  * F(fi.i)  at  r — 0 (9km) 

where  G'(»)  is  the  radiative  equilibrium  condition  and  £(«)  is  the  iropopause  geopolential 


height  field  derived  from  climatological  data  as  discussed  in  Paper  I.  The  tropopause  forcing  of 
planetary  waves,  FiH.i),  is  identical  to  that  used  by  Matsuno  (1971)  and  is  applied  only  to  the 
winter  hemisphere. 
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III.  Nunirrical  Results  for  In  damped  Warmings  in  an 
Isothermal  Atmosphere 

a.  m = I 

Our  first  numerical  experimeni  was  an  aitempi  to  reproduce  Maisuno  s (1971)  results 
An  isothermal  atmosphere  of  256K  was  chosen  to  approximate  the  static  stability  of  the  late 
winter,  high  latitude  stratosphere  The  initial  wind  held  (I'lg  la)  was  similar  to  Maisuno  s and 
damping  was  set  to  zero  to  allow  only  eddy  stress  lo  alter  the  wind  contiguration  In  the 
mechanistic  model  used  by  Matsuno  the  generation  of  a SSW  event  was  accomplished  by  the 
switch-on  of  a geopoieniial  height  perturbation  over  a 5 day  period  to  an  asymptotic  value  ol 
.700  gpm  at  the  lower  boundary  of  300  mb.  With  the  same  lower  boundary  condition,  we 
obtained  the  results  for  am-  1 warming  shown  as  time-height  contours  lor  <*  , i7,  T,  and  T in 
Fig  2 The  associated  mean  zonal  wind  profiles  with  latitude  at  10  and  25  days  into  the  simula- 
tion are  presented  in  Fig  3 

From  Fig  2c  it  is  clear  that  a strong  double  warming  (AT  ~45K)  is  produced  after  15-2.' 
days  of  integration  at  30  and  60  km  and  is  accompanied  by  strong  cooling  in  the  upper  meso- 
sphere The  zonal  wind  field  develops  two  easterly  regions  independently  at  100km  and  30km 
associated  with  the  double  warming  which  expand  downward  and  upward  to  pinch  off  the  polar 
night  let  (cf  Figs  2b  and  3b).  The  general  features  of  this  m = 1 double  warming  resemble 
closeh  Maisuno’s  double  warming  and  would  tend  to  confirm  his  calculations. 

However  the  upper  level  warming  does  not  occur  as  a result  of  critical  level  descent  from 
the  mesosphere  as  Maisuno  hypothesized.  An  examination  of  Fig.  2a  reveals  strong  transient 
oscillations  in  the  middle  atmosphere  The  switch-on  of  stationary  wave  forcing  at  the  lower 
boundary  launches  a spectrum  of  traveling  planetary  waves  that  propagate  in  both  longitudinal 
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directions.  The  eastward  moving  waves  are  probably  unimportant  since  they  are  evanescent  in 
regions  where  the  mean  zonal  wind  is  seen  as  easterly  and  decay  when  the  lower  boundary  fore- 
seen as  easterly  and  decay  when  the  lower  boundary  forcing  becomes  steady.  However,  some 
ing  becomes  steady.  However,  some  of  the  westward  moving  waves  can  become  trapped 
between  the  polar  night  jet  and  the  ground  (or  the  lower  boundary  in  our  case).  In  particular 
Clark  and  Schoeberl  (1979)  found  a strong  resonantly  trapped  planetary  wave  with  a period  of 
IS  days  for  m = 1.  These  traveling  planetary  waves  will  move  in  and  out  of  phase  with  the  sta- 
tionary planetary  wave  to  produce  the  oscillation  in  wave  amplitude  displayed  in  Fig.  2a. 


To  illustrate  this  point  more  rigorously  let  us  examine  simplified  versions  of  Eqs.  (5)  and 
(6),  respectively,  on  the  beta  plane 
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where  g represents  potential  vorticity  (Dickinson,  1969;  Holton,  1975).  The  wave  motion  can 
be  separated  into  a stationary  component  (s)  with  exp  Um\)  dependence  and  a traveling  com- 
ponent (t)  with  exp  (im\  + kti)  dependence  where  K is  longitude.  A straightforward  calcula- 
tion with  Eqs.  (8)  and  (9)  yields 
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where  v,  , denotes  the  amplitude  and  phase  of  the  northward  wave  velocity.  If  v,  and  v,  are  in 
phase,  then  Eq.  (10)  reduces  to 
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Thus  the  mean  zonal  wind  and  temperature  fields  forced  by  stationary  wave-traveling  wave 
i.'tteraclion  oscillate  at  the  frequency  of  the  traveling  wave.  Although  the  mean  zonal  wind  is 
both  accelerated  and  decelerated  during  each  oscillation,  the  time  averaged  acceleration  is  zero 
except  at  a critical  level  in  agreement  with  wave-zonal  fiow  interaction  theorems. 
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Some  of  the  results  of  the  numerical  experiment  presented  m Fig  2 can  now  be  clearly 
interpreted  in  terms  of  transient  w-ave-staiionury  wave  interaction  After  the  upward  propagai- 
mg  stationary  wave  front  initially  decelerates  the  zonal  wind  as  predicted  by  Matsuno  (IdTII 
and  r.yu  (1974),  the  deceleration  is  continued  m the  lower  stratosphere  by  the  stationary 
wave-traveling  wave  interaction  expressed  by  Flq.  (10),  however  once  a critical  level  is  esta- 
blished for  any  wave  component,  further  mean  zonal  wind  deceleration  occurs  at  the  wave 
absorbing  critical  level,  independently  of  wave-wave  interaction  A detailed  comparison  of  I'lgs 
2a,  2b,  and  2c  reveals  a close  correlation  of  the  wave  amplitude  with  u and  /'with  a time  delay 
consistent  with  F'q.  (II)  For  example,  note  the  rapid  change  from  strong  easterlies  to  weak 
westerlies  at  b5km  during  the  period  days  2,7-28  that  follows  the  rapid  decrease  m maximum 
wave  amplitude  at  the  same  location  for  the  period  days  20-25.  It  should  be  pointed  out  that 
two  transient  waves  with  different  frequencies  could  interact  in  an  analogous  manner  to  the  sta- 
tionary wave-traveling  wave  interaction.  A precise  determination  of  the  amplitudes  and  fre- 
quencies of  the  transient  waves  cannot  be  made  because  the  mean  zonal  How  vacillation  con- 
tinuously shifts  the  fundamental  frequencies  of  these  waves  and  generates  additional  com- 
ponents. 

While  the  upper  level  warming  at  bOkm  in  Fig  2c  more  clearly  demonstrates  transient 
wave  effects  in  addition  to  critical  level  absorption  of  the  stationary  wave,  the  lower  stratos- 
pheric warming  appears  to  be  controlled  by  critical  level  deceleration  once  easterlies  have  been 
established  (cf.  Fig.  2b).  Only  a very  small  amplitude  transient  disturbance  is  superimposed  on 
this  warming.  During  the  evolution  of  this  warming  the  easterly  wind  line  descends  approxi- 
mately 15  km  consistent  with  Matsuno’s  (1971)  explanation  ol  critical  level  interaction. 


Ilirota  (1968)  and  .Schoeberl  (1978)  suggested  that  thermal  planetary  wave  may  be  an 
tmportant  observational  component  of  the  .S,SW  structure  In  Fig  2tl  the  amplitude  ol  / at 


NKl.  Ml  MORANIH'M  Rl  PdRT  4<)IJ 


60°N  is  shown  to  be  comparable  to  variations  in  T at  the  polar  stralopause.  In  the  lower  level 
warming  its  amplitude  is  only  one-fourth  of  the  T variation. 

b.  m “=  2 

The  development  of  the  m = 2 warming  is  strikingly  different  from  the  m = 1 warming 
as  can  be  seen  from  Figs.  4 and  5.  For  m = 2 the  warming  develops  explosively  after  20  days 
with  a 60K  temperature  increase  in  4 days  at  35  km  to  a maximum  of  70K,  almost  as  large  as 
that  obtained  by  Matsuno  (1971).  We  also  obtain  a relatively  weak  secondary  warming  at 
60km;  otherwise,  the  basic  features  of  our  m = 2 warming  are  very  similar  to  Matsuno’s 
(1971)  results.  The  m = 2 wave  amplitude  (Fig.  4a)  is  much  weaker  than  the  m — 1 wave 
amplitude  (Fig.  2a)  during  the  first  15  days.  But  from  day  15  the  wave  amplitude  increases 
dramatically  in  the  upper  stratosphere  and  reaches  a maximum  at  day  21,  followed  by  a rapid 
increase  in  wave  amplitude  in  the  upper  mesosphere  to  a comparable  maximum  on  day  24. 
The  zonal  wind  deceleration  is  extremely  rapid  after  day  17.  and  approximately  day  20  easterlies 
appear  nearly  simultaneously  throughout  the  upper  stratosphere  and  mesosphere  (Fig.  4b). 

From  Fig.  5 the  behavior  of  the  m = 2 warming  can  be  understood.  In  contrast  to  a 
m - 1 warming  which  generates  an  equatorward  moving  critical  level  (Fig.  3b),  the  easterly 
wind  region  in  the  m =-  2 event  advances  poleward  from  the  tropics  as  shown  in  Fig.  5b.  The 
warming  occurs  suddenly  as  the  poleward  moving  critical  line  in  the  stratosphere  intercepts 
sufficient  vertically  propagating  wave  energy  to  produce  easterlies  over  most  of  the  middle 
atmosphere  (Fig.  5c).  Figs.  5b  and  5c  also  reveal  a separate  descending  critical  level  from  the 
equatorial  mesosphere  which  is  formed  by  the  transient  wave  pulse  launched  by  the  switch-on 
lower  boundary  condition  and  contributes  to  the  upper  level  warming  (Fig.  4c).  On  the  basis  of 
other  experiments,  the  development  of  the  equatorial  easterly  region  appears  to  be  "uiependent 
of  the  initial  wind  configuration  in  the  tropics  and  is  a consequence  of  the  divergence  of  eddy 
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momentum  (lux  (i/v  ) in  equatorial  regions.  Since  v'  <r  m<t>\  this  flux  is  twice  as  large  for  m - 
2 as  for  w - 1.  when  <t>'  is  comparable  These  results  are  consistent  with  Holton’s  (1976)  prim- 
itive equation  simulation  of  a m - 2 warming  which  likewise  produced  poleward  advancing  crit- 
ical levels  from  equatorial  regions. 


The  postwarming  period  lor  ni  — 2 does  not  display  the  pronounced  oscillations  which  are 
characteristic  of  traveling  planetary  waves.  I he  lack  o(  ni  "■  2 traveling  waves  or  free  modes  is 
consistent  with  the  recent  calculations  of  Clark  and  Schocberl  (1979),  who  found  no  strong  m 

- 2 resonances  for  westward  propagating  waves  with  periods  between  5 and  50  days.  Unlike  m 

- 1 waves,  m - 2 planetary  waves  experience  greater  vertical  trapping  by  stratospheric  winds 
(Matsuno,  1970;  Schoeberl  and  Geller,  1976). 


In  Fig.  4d  some  characteristics  of  the  m « 1 thermal  planetary  wave  are  illustrated. 
Although  its  amplitude  is  substantial  (—  20K).  the  large  variation  of  f would  probably  swamp 
its  signature  in  noisy  observational  data.  Unlike  the  m - 2 thermal  planetary  wave,  the  m - 2 

thermal  wave  has  large  amplitude  spread  over  the  30-70  km  region.  Since  T <r  the  ther- 

dz 

mal  wave  structure  for  ni  - 1 and  2 warming  is  illustrative  of  vertical  trapping  of  wave  energy 
by  the  polar  night  jet.  Greater  trapping  generates  more  vertical  variation  in  4,'  and  hence  larger 
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IV.  NUMERICAL  RESULTS  WITH  DAMPING 

While  the  results  in  the  previous  section  can  illuminate  certain  dynamical  processes,  they 
do  not  apply  to  the  real  atmosphere  where  thermal  wave  energy  may  be  radiated  to  space  and 
wave-wave  interaction  can  ultimately  lead  to  dissipation  of  wave  momentum  Similarly  m the 
absence  of  any  wave  forcing  the  /onally  averaged  circulation  would  return  to  the  prewarming 
configuration  by  these  dissipative  processes,  l-'or  lack  of  a precise  theoretical,  yet  computation- 
ally economical  description  of  these  processes,  we  adopt  the  strategy  of  Paper  I and  represent 
these  dissipative  processes  by  Newtonian  cooling  and  Rayleigh  friction  with  the  same  numerical 
values  for  both  the  waves  and  mean  How  as  given  in  Paper  I It  should  be  remembered  that  tbe 
height-dependent  damping  coellicients  determine  the  relaxation  time  of  the  mean  How  after 
cessation  of  wave  eddy  stresses. 

In  this  section  integrations  for  m - I and  rn  - 2 warmings  are  reported  which  were  per- 
formed with  damping  and  the  initial  wind  field  in  Fig  lb  (cf  Section  ID  Table  1 contains  a 
summary  of  all  the  integration  reported  in  this  paper.  The  post  warming  phase  of  these  integra- 
tions will  be  discussed  in  Section  VI. 

a.  m — 1 

The  /n  — 1 results  are  displayed  in  Figs.  6 and  7 and  should  be  compared  with  Figs.  2 and 
3 to  determine  the  impact  of  damping.  Figure  6a  clearly  shows  that  damping  reduces  the  wave 
amplitude  by  approximately  a factor  of  2 and  smooths  out  the  wave  oscillation  to  a regular 
period  of  14  days  (Figs.  2a  and  6a).  The  period  of  this  traveling  wave  which  is  resonantly 
trapped  between  the  polar  night  jet  and  the  lower  boundary  is  sensitive  to  the  zonal  wind 
configuration  and  the  boundary  condition  (Clark  and  Schoeberl,  1979).  A comparison  of  Figs. 
6b  and  2b  shows  that  there  is  less  distortion  of  the  mean  flow  by  wave  eddy  stresses  when 
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damping  is  included  Thus  the  period  ol  the  traveling  wave  should  be  less  variable  as  our 
numerical  results  indicate 

With  the  presence  o(  damping  only  very  weak  easterlies  are  formed  at  —100  km  (the 
“mesospheric"  critical  level)  from  the  transient  wave  pulse  Tliesc  easterlies  cannot  be  sus- 
tained by  the  stationary  planetary  wave  stresses  along  the  critical  level  against  the  restoring 
force  of  damping  which  is  large  near  the  upper  boundary.  In  addition  very  little  wave  ampli- 
tude reaches  this  region  after  day  .15  as  a consequence  of  attenuation  as  it  propagates  through 
the  stratosphere  The  stratospheric  critical  level  at  .15  km  does  not  develop  until  — 50  days  and 
then  only  for  a short  duration 

■M  higher  altitudes,  small  warmings  initially  develop  at  — 65  km  (Fig  6c)  clearly  not 
associated  with  any  critical  level  but  exhibiting  a close  correlation  with  oscillations  in  wave 
amplitude  (Figs  6a  and  6b)  and  disappears  after  .10  days  as  the  traveling  waves  decay  The 
major  warming  at  .15  km  ts  still  in  its  developmental  pha.se  at  this  time  land  does  not  reach  its 
peak  of  — 22  K until  day  551  With  damping  included  this  warming  lakes  twice  the  time  and  is 
a factor  of  2 weaker 

A close  examination  of  the  numerical  results  indicates  that  this  warming  is  not  due  to  crit- 
ical level  elTecis  or  stationary  wave-transient  wave  interaction  but  rather  damping  With  damp- 
ing present  wave-mean  flow  interaction  can  lead  to  a net  deceleration  of  the  mean  zonal  wind 
f-'or  the  linear  damping  used  in  our  model  deceleration  by  stationary  waves  is  strongest  where 
the  wave  amplitude  is  large  and  the  zonal  wind  is  weak  (cf.  Kq.  (5))  The  warming  occurs  in  a 
region  where  the  restoring  force  of  local  damping  on  the  mean  flow  is  not  sufficient  to  over- 
come the  planetary  wave  stresses  produced  by  the  convergence  of  wave  energy  and  momentum 
Fventually  the  zonal  flow  is  decelerated  to  easteiiies. 
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The  thermal  planetary  wave  reaches  a maximum  of  ~ 20  K m the  region  of  the  small, 
high  altitude  warming  produced  by  stationary  wave-traveling  wave  interaction  and  is  primarily  a 
traveling  thermal  disturbance  Strong  damping  of  these  traveling  oscillations  m the  low  alti- 
tude, major  warming  region  prevents  the  growth  of  large  amplitude  thermal  planetary  waves  in 
this  type  of  warming  event 

b 01-2 

The  introduction  of  damping  has  the  same  eH'ecis  on  a oi  — 2 warming  as  discussed  above 
for  the  m - I warming.  The  warming  lakes  approximately  twice  as  long  to  develop  as  without 
damping  and  its  amplitude  is  almost  decreased  by  a factor  of  2 (Figs  8 and  9).  Damping  also 
slows  the  explosive  growth  of  the  m - 2 warming.  The  development  of  the  poleward  advanc- 
ing critical  level  initially  formed  in  the  eguatorial  stratosphere  is  considerably  retarded  as  damp- 
ing acts  to  restore  the  original  westerly  wind  configuration  (Fig.  9).  The  eddy  stresses  must  stis- 
taut  the  easterlies  as  well  as  advance  the  critical  level  poleward.  The  equatorial  mesospheric 
critcal  level  never  really  develops  or  advances  (Figs.  9b-9d),  since  upward  propagating  waves 
are  blocked  by  the  stratospheric  critical  level  and  damping  is  strong  enough  at  these  height;,  to 
restore  westerlies. 

As  might  be  expected  from  the  first  experimeni,  neither  the  wave  amplitude  (Fig,  8a)  nor 
the  thermal  wave  (Fig.  8d)  exhibit  the  periodic  oscillations  so  evident  in  m - 1 warmings.  In 
addition,  the  thermal  wave  amplitude  (T  ) is  considerably  less  than  the  increase  of  f at  the 
peak  of  the  warming.  But  note  that  prior  to  day  .15,  T is  comparable  to  the  change  in  T. 
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V.  OTHER  NUMERICAL  RESULTS 

In  addition  to  the  introduction  of  damping  there  are  other  input  parameters  that  can  play 
an  important  role  in  the  evolution  of  warming  From  our  previous  discussion  the  strength  of 
the  lower  stratospheric  westerly  winds  is  a strong  candidate.  Since  an  induced  secondary  circu- 
lation is  an  important  component  in  the  stratospheric  response  to  eddy  heal  transport  we  also 
anticipate  that  the  static  stability  ol’  the  stratosphere  could  potentially  be  Important. 

a.  Nonisothermal  Basic  State 

Numerical  experiments  were  performed  for  m — 1 and  m - 2 warmings  with  the  static 
stability  evaluated  from  the  globally  averaged  CIRA  (1972)  temperature  profile.  The  results  are 
given  in  Table  1.  The  overall  effect  was  minor,  but  the  major  stratospheric  warming  occurred 
4-5  days  earlier  for  each  zonal  harmonic.  In  the  case  of  m - 2,  the  warming  was  8 K more 
intense.  Below  28km  the  static  stability  of  the  nonisothermal  atmosphere  was  less  than  that  of 
isothermal  atmosphere.  A more  vigorous  induced  secondary  mean  circulation  was  generated 
from  the  greater  available  potential  energy  in  the  mean  flow.  The  net  effect  was  a stronger 
warming  as  measured  by  Af  and  a faster  evolution  of  the  warming. 

b.  Lower  Stratospheric  Winds 

To  investigate  the  role  of  lower  stratospheric  winds  an  isothermal  atmosphere  was  used 
and  M was  set  equal  to  zero  in  Eqs.  (1)  and  (6).  This  decreased  the  initial  wind  configuration  in 
the  lower  stratosphere  from  the  strength  in  Fig.  lb  to  the  weak  westerlies  illustrated  in  Fig  Ic. 
From  the  summary  in  Table  1.  we  note  that  the  polar  temperature  rise  (f)  for  both  w - 1 and 
m - 2 is  reduced  by  a factor  of  three  due  to  weaker  winds.  This  result  is  consistent  with  the 
theoretical  analysis  of  Dickinson  (1969b)  who  showed  that  the  transmission  of  planetary  waves 
in  the  presence  of  strong  damping  and  weak  westerlies  is  very  low.  The  increased  absorption  of 
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wave  energy  near  the  lower  boundary  does  not  initiate  a warming  in  that  region  because  of  the 
imposed  6 boundary  conditions  at  300  mb  For  the  m - 2 simulation  the  equatorial  easterly 
regions  advances  poleward,  as  belore,  and  establishes  an  easterly  regime  throughout  the  middle 
stratosphere  However  very  little  wave  energy  penetrates  through  the  weak  westerlies  to  pro- 
vide a substantial  circulation  reversal.  Also  the  critical  level  cannot  descend  to  a lower  altitude 
tor  reasons  discussed  above  and  absorb  more  wave  energy.  For  m - 1 the  severe  attenuation 
of  planetary  wave  energy  inhibits  the  development  of  the  pronounced  periodic  oscillations  in 
the  mesosphere  which  were  characteristic  of  previous  m — 1 simulations.  For  example  at  70 
km  12b  dm  at  day  15,  but  only  5 dm  for  its  next  peak  on  day  30. 

c fri  - 3 warming 

An  attempt  was  made  to  generate  a m - 3 warming  even  though  ni  - 3 warmings  have 
never  been  observed  After  bO  days  of  integration  no  major  warming  was  produced  and  the 
results  did  not  suggest  further  integration  time  was  warranted  Planetary  waves  with  high  zonal 
wavenumbers  (w  > 2)  are  strongly  trapped  by  the  polar  night  let  and  become  increasingly 
evanescent  throughout  the  bulk  of  the  stratosphere  (Charney  and  Prazin,  l*7bl)  Thus  the  ina- 
bility to  simulate  a iii  - 3 warming  is  not  surprising. 

d Mesospheric  cooling 

For  both  m - 1 and  m — 2,  strong  mesospheric  cooling  accompanies  the  strong  stratos- 
pheric warming  as  a result  of  the  induced  secondary  circulation  and  associated  adiabatic  cooling 
in  response  to  eddy  heat  transport  to  the  polar  stratosphere  This  effect  is  clearly  seen  in  Fig- 
ures 6c  and  8c  For  w - 1 the  maximum  cooling  generated  near  the  pole  is  12  K in  the  80-100 
km  region  while  for  m - 2.  the  me,sosphere  cools  22  K in  the  same  region  These  values  are 
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in  reasonably  good  agreement  with  those  of  Labitzkes  (1972)  who  observed  cooling  above  the 
stratopuase  on  the  orders  of  20  K at  high  latitudes  during  a major  warming. 
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VI.  LONG  TIME  INTEGRATIONS;  THE  ASYMPTOTIC  STATE 

The  previous  studies  on  straiosphene  vacillalion  eseles  by  Holton  and  Mass  (1976)  and 
Paper  1 led  us  to  carry  out  miegraiions  for  m - 1 and  ni  - 2 to  120  and  180  days,  respectively, 
vciih  asymptotic  forcing  of  300  gpm  at  300  mb  Figures  10  and  II  present  the  representative 
values  of  <t>  and  i7  at  120  days  for  m - 1 and  140  days  for  rn  ^ 2 At  those  times  m - 1 and 
fn  “ 2 are  in  post  warming  periods  The  structure  ot  ft>  for  m **  I and  m » 2 show  similar 
maxima  in  the  mesosphere  and  secondary  maxima  in  the  lower  stratosphere  (Fig  10) 

A comparison  of  Fig  II  with  the  initial  wind  configuration  in  Fig  lb  shows  that  the  polar 
night  jet  strength  is  reduced  by  approximately  8-9  ms  ' at  its  maximum  through  the  action  of 
planetary  wave  stresses  However  the  largest  changes  occur  in  the  lower  stratospheric  westerlies 
which  are  decelerated  by  — 20  ms  ' through  eddy  heat  transport  by  planetary  waves  into  the 
polar  stratosphere  For  m — 1 the  polar  temperature  increases  by  20K  at  30  km,  but  declines 
by  8K  at  80  km  as  a result  of  the  adiabatic  cooling  associated  with  the  induced  secondary  meri- 
dional circulation  Similarly  for  m - 2 planetary  waves  the  polar  temperature  increases  by  15K 
at  30  km,  and  decreases  by  9K  at  80  km  Thus  planetary  waves  tend  to  cool  the  polar  meso- 
sphere rather  than  heat  it  as  Cireen  (1972)  suggested  The  computed  effects  of  these  waves  are 
somewhat  larger  than  might  be  expected,  since  we  have  used  300  m forcing,  about  twice  the 
observed  monthly  mean  wave  amplitude  at  300  mb  for  January  (van  Loon  et.  al  . 1973) 

In  Paper  1 steady  state  solutions  for  planetary  waves  iiiic  ’ -g  with  the  mean  zonal  flow 
were  found  for  the  m - 2 wave  For  the  m - I wave  a mesospheric  critical  level  was  formed 
after  the  first  iteration  and  rendered  the  iterative  procedure  divergent  in  our  search  for  a steady 
state  solution  Our  present  calculations  indicate  that  a steady  state  m - 1 solution  does  exist, 
but  can  only  be  obtained  by  time  dependent  calculations 
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llolion  and  Mass  (1*>76)  and  Holton  and  Dunkerton  (1978)  tiave  noted  that  the  stratos- 
pherit  vatdiaiion  c>t'le  ol  planetarv  »ase  amplitude  and  mean  circulation  can  occur  even  il  the 
»ase  lorcinR  at  the  lo*er  boundary  is  steads  The  appearance  of  a second  ni  - 2 warming  after 
lias  I It)  and  a third  warming  at  da>  Ih8  during  the  integration  is  certainly  suggestive  of  a vacil- 
lation i.vcle  I rom  I ig  12  note  that  each  succeeding  in  - 2 warming  is  slightly  weaker  than  the 
previous  one  Mthough  this  integration  was  extended  to  b months  it  would  appear  that  an 
ec|uilihrium  solution  is  aitainahle  onlv  alter  more  than  a year's  integration  The  adjustment  pro- 
cess in  mulipic  "I  2 warmings  is  initimatelv  coupled  to  the  lateral  motion  of  the  critical  level 
Kut  this  motion  was  not  alh>wed  bv  the  Holton  and  Mass  (197b)  model  Thus  tt  is  dtfticult  to 
compare  our  model  with  its  much  moie  complex  treatment  of  the  mean  zonal  flow  with  their 
results 

1 nlike  Holton  and  Dunkerton  (!9’8)  who  suggest  the  mechanism  for  the  vacillation  cycle 
IS  transience  from  the  initial  switch-on  forcing,  we  find  it  more  plausible  to  attribute  the  vacilla- 
tions to  a cyclic  formation  of  critical  levels,  followed  hy  restoration  of  the  westerlies  through 
damping,  and  a suseguent  increase  in  planetarv  wave  transmission  and  eddy  heat  transport  to 
re-establish  another  critical  level  After  each  critical  level  is  formed  m the  equatorial  regions  it 
moves  poleward  until  it  intercepts  all  of  the  upward  propagating  wave  energy  In  our  model  the 
critical  level  cannot  descend  to  the  lower  boundary  for  reasons  discussed  in  Section  \’.  How- 
ever, in  the  real  atmosphere  the  critical  level  could  descend  into  the  troposphere  to  a level 
which  would  cut  off  the  source  of  planetarv  wave  forcing.  Thus  our  artilkial  lower  boundary  at 
the  tropopause  is  not  unrealistic  Once  the  easterly  wind  regime  is  established,  the  topside  of 
the  easterly  wind  region  begins  to  erode  as  the  mean  How  is  damped  and  the  easterly  region 
eventually  becomes  thin  enough  to  permit  penetration  of  planetarv  waves  into  the  upper  strato- 
sphere albeit  severely  attenuated  Damping  continues  to  slowly  restore  the  winds  to  the 
prewarmmg  configuration,  at  which  point  the  entire  process  repeats  itself 
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Vll.  CONCLLDING  REMARKS 

The  linear  damping  used  m our  calculations  should  be  accurate  in  the  initial  development 
of  a warming  However,  once  the  planetary  waves  and  the  changes  in  the  zonally-averaged  cir- 
culation reach  finite  amplitude,  this  approximate  representation  of  damping  probably  breaks 
down  and  leads  to  unrealistic  estimates  of  the  actual  effects  of  radiative  damping  and  frictional 
dissipation  In  the  case  of  Rayleigh  friction  there  may  also  be  a contribution  of  unknown  mag- 
nitude that  represents  planetary  wave  stresses  on  the  mean  flow  and  is  inappropriate  for  damp- 
ing of  planetary  waves.  Thus  the  numerical  simulations  that  were  performed  with  and  without 
damping  should  be  regarded  as  extreme  bounds  on  the  actual  development  of  SSW  events  in  a 
realistic  atmosphere.  In  spite  of  this  limitation  our  model  clearly  delineates  many  essential 
dynamical  processes  that  control  and  differentiate  m = 1 and  m — 1 warmings. 

The  development  of  a stratospheric  warming  in  the  real  atmosphere  depends  critically  on 
two  factors;  (1)  the  strength  of  the  westerly  winds  in  the  lower  stratosphere  and  (2)  the  magni- 
tude of  wave  damping  in  the  same  region.  These  two  factors  govern  the  transmission  of  plane- 
tary wave  energy  to  the  less  dense  regions  of  the  stratosphere  and  mesosphere  where  warmings 
originate.  Large  transmission  requires  wind  velocities  greater  than  25  ms*'  everywhere  at 
midlalitude  or  negligible  damping  (~  100  day*').  Thus  we  expect  that  major  warmings  will 
develop  only  when  the  prewarming  lower  stratospheric  winds  are  large.  An  observational  check 
of  this  key  theoretical  result  is  difficult  because  most  stratospheric  wind  data  is  monthly  aver- 
aged. Over  this  time  period,  the  possibility  that  prewarming,  warming,  and  postwarming  data 
remain  distinct  is  remote. 

Damping  has  two  obvious,  important  effects  on  the  evolution  of  both  m = 1 and  in  — 2 
svarmings.  It  can  suppress  the  maximum  amplitude  that  a warming  attains  and  slow  its  develop- 
ment. These  effects  are  consistent  with  Matsuno's  (1971)  and  Geisler’s  (1974)  earlier  studies 
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(■Dr  a tti  “ 1 Ncarniiim;  damping  can  be  the  most  imporiani  factor  that  results  m the  convergence 
of  planetars  ssase  heat  and  momentum  llu.xes  to  decelerate  the  zonal  wind,  heat  the 
stratosphere,  and  precipitate  the  warming  Hut  the  most  distinctive  feature  of  a tti  = 1 warming 
IS  pronounced  oscillations  with  a period  of  approximately  two  weeks.  These  oscillations  are  the 
result  of  resonantK  trapped  westward  propagating  planetary  waves  (Clark  and  Schoeberl,  1979) 
moving  in  and  out  of  phase  with  the  tropospherically  forced  stationary  planetary  wave  .-ks  a 
result  of  the  interaction  of  the  stationary  and  traveling  waves,  the  mean  zonal  How  oscillates  at 
the  traveling  wave  frequency  This  oscillation  can  reach  sufficient  amplitude  at  times  to 
decelerate  the  zonal  flow  during  a cycle  to  easterlies  and  thus  create  a critical  level  In 
Matsuno's  (1971)  model,  the  upper  part  of  the  double  warming  is  generated  by  these  oscilla- 
tions However,  even  if  a critical  level  is  not  formed,  the  temperature  increase  associated  with 
the  oscillation  has  the  characteristic  magnitude  of  a warming  and  may  be  regarded  as  the  initial 
phase  of  a m = 1 warming  Damping  plays  the  most  important  role  in  the  later  phase  of  a 
m =■  1 warming  and  can  lead  to  a stronger  temperature  increase  if)  than  the  initial,  oscillatory- 
phase 

Although  formation  of  a mesospheric  critical  level  is  not  a required  feature  of  a warming, 
the  development  and  propagation  of  critical  levels  is  central  to  the  evolution  of  the  sudden 
warmings  Initially  in  a w ” 1 event  the  critical  level  forms  in  the  polar  mesopause  region  and 
advances  equatorward  and  downward  in  its  development.  In  sharp  contrast  the  mesospheric 
critical  level  in  a w - 2 warming  initially  develops  m the  equatorial  mesopause  region  and 
advances  poleward  This  difference  is  likely  due  to  the  greater  refraction  of  the  wave  energy- 
flow  toward  the  equator  by  the  mean  wind  field  for  m = 2 than  for  m = 1 (Matsuno,  1970). 
Similar  remarks  about  the  propagation  of  critical  levels  applies  to  their  development  in  the  stra- 
tosphere Thi  »Fi  - I critical  level  initially-  form  in  polar  regions  and  grows  equatorward, 
whereas  w — 2 critical  level  originates  at  equatorial  latitudes  and  progresses  poleward.  Since  the 
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m - I planelary  wave  ean  propagate  with  larger  amplitude  into  high  latitude  regions  of  the  stra- 
tosphere than  a in  - 2 wave,  the  convergence  of  eddy  heat  llux  lor  in  — 1 occurs  in  the  more 
contined  polar  region  and  leads  to  a greater  local  heating  rate  than  the  midlatiiude  conlined, 
III  —2  wave  In  addition,  the  greater  trapping  of  the  in  — 2 wave  relative  to  the  in  — I wave  by 
the  jei  also  reduces  the  eddy  heat  llux  to  high  latitudes.  Thus  the  initial  critical  level  formation 
occurs  in  the  polar  regions  for  w - I.  Near  the  equatorial  regions  both  in  - I and  in  - 2 
decelerate  the  mean  zonal  How  as  a result  of  momentum  llux  divergence  but  in  - 2’s  ell'ect  is 
substantially  l.irger  since  v'  is  twice  as  large  for  similar  amplitude  waves  Thus  the  in  — 2 
warming  is  dominated  by  stronger  equatorial  deceleration  of  the  Ilow  and  a polewaril  moving 
critical  layer. 

I'he  other  distinctive  feature  of  a in  ~ 2 warming  (in  addition  to  its  poleward  advancing, 
critical  level)  is  the  sudden,  explosive  character  of  its  development  even  in  the  presence  of 
damping  Haseil  on  our  arguments  that  damping  at  finite  amplitudes  is  probably  unrealistic  in 
our  calculations,  a growth  rale  intermediate  between  the  rates  exhibited  in  I'lgs  4c  and  Sc  may 
be  appropriate.  I'rom  our  calculations  we  would  infer  that  in  - 2 planelary  waxes  must  play  an 
essential  role  in  observed  warming  events  that  develop  suddenly  l.abil/ke  (1477)  noted  that  a 
in  - 2 planelary  wave  appears  to  be  ineapable  of  generating  a ma.ior  warming  by  itself,  but 
Sehoeberl  (1478)  points  out  that  the  interaction  of  both  in  - 1 and  in  - 2 is  observed  m most 
observed  warmings  I hat  is,  in  — 2 doc.s  play  some  role  in  warming  development,  even  though 
the  amplitude  of  in  - I is  strongest  just  prior  to  the  warming. 

C'onsisleni  with  observations,  our  numerieal  simulations  also  prediei  realisiie  values  of 
polar  mesospheric'  eooling  eoineideni  with  a polar  siraiospherie  warming.  This  is  a direei  conse- 
quenee  of  the  indueed  secondary  circulation  in  response  to  eddy  heal  transport  to  the  polar  stra- 
tosphere which  results  in  adiabatic  cooling  in  the  mesosphere  l.abiizke  (1472)  and  lliroia  and 
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Hau\oU  U'477i  h(»\,o  (.icvluccil  nicsosphcriv  cooling  in  pol.ii  regions  Ironi  aiiiilysis  of  wanning 
ocenis  Based  on  our  long  lime  iniegiaiion  siudies  we  find  dial  planeiaiv  waves  ihrough  ilie 
induced  secondary  circulaiion  cause  a net  cooling  of  the  polar  mesosphere 

1 rom  an  obsei vaiional  point  of  view  our  calculations  also  predict  a variety  of  distinctive 
features  that  have  vet  to  he  confirmed  with  observational  data  1 hev  include  the  need  for 
stioiig  westerlies  m the  lower  stratosphere,  the  dilVerent  latitudinal  regions  of  ni  - I and  in  - 2 
critical  level  formation,  the  pronounced  in  — 1 oscillations  with  a period  - 15  days  (Madilen, 
l‘)7s,  ('lark,  .iiul  Schoeberl,  Id?')),  the  strong  thermal  planetary  wave  imbedded  m in  — 1 
warmings  (Kirota,  IdbX.  I abit/ke.  Id??,  Schoeberl,  ld78l,  and  periodic  cntical  level  foimalion 
111  in  *■  2 warmings 

f inally  our  results  suggest  that  sudden  stratospheric  warmings,  are  a necessarv  part  of  the 
overall  adivistment  process  wbeiehy  the  atmosphere  responds  to  the  strong  radiative  cooling  m 
polar  regions  during  the  long  winter  night  Normally  the  synoptic  scale  eddies  ihrough  well 
known  barocimie  processes  in  the  troposphere  and  lower  stratosphere  transport  heat  polewards 
ti'  counter  balance  radiative  cooling  ,\s  winter  progresses  the  polar  siratospheie  cools  dramali- 
callv  and  forms  the  intense  polar  night  lei  to  maintain  thermal  wiml  balance  .-ks  this  vortex 
evp.iitvls  downward  and  mtensilies,  the  svnoptic  scale,  barocimie  waves  eventuallv  become 
ev.mesceiil  in  the  lower  stratosphere  and  no  longer  transport  heat  poleward  to  maintain  polar 
sir.'tosphenc  iemper,iiures  Hiiwever,  m the  presence  ol  a strong  polar  night  let  the  large  scale 
p!  inetarv  waves  can  penetrate  through  the  lower  stratosphere  with  large  transmission  m spite  of 
damping  .ind  transport  heal  (loleward  fhus  the  burden  of  maintaining  pidar  stratospheric  tem- 
per.iiures  ,ig,iinsi  radiative  cooling  shifts  from  the  synoptic  sciile  to  planetary  scale  waves  The 
sudden  warming  is  ,i  product  ol  this  shift 
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\pp«ndix  A — Notation 


In 


P 


, where  p 


pressure.  p„  - 300  mb 


eastward  and  northward  wind  velocities,  respectively 

llL 

ill 

earth’s  radius 

angular  rotation  rate  of  the  earth 

dry  air  gas  constant 

latitude 

zonal  wave  number 
longitude  (measured  east) 
geopntential  height 

zonally  averaged  heating  rate  (departure  from  global  average) 
Rayleigh  friction  coefficient 
Newtonian  cooling  coefficient 


+ — ; . where  <T>  is  the  globally  averaged  temperature  and  c„ 

i;  bz 

specific  heat  at  constant  pressure 
external  momentum  sources  and  sinks 
{kjlil)' 

globally  averaged  variable 


zonally  averaged  variable  - < > 
perturbation  from  ( ) 


3^  _1_ 

bS 

■ ^ 1 

bS 

b'S 

4 

bz 

^Ts 

bz 
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Appendix  B — Numerical  Model 


Defining 


ti/  - <■  - S ■ <t> 


we  can  transform  Eqs.  (5)  and  (6)  to 
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Kquations  (Bl)  and  (B2)  are  a coupled  system  of  lime-dependent  elliptic  partial 
differential  equations  which  may  be  solved  by  standard  numerical  techniques.  Both  ^ and 
L\(4>)  + Z.  1 „(i/( ) are  required  to  compute  the  non-linear  forcing  term  in  Eq.  (Bl ).  An  explicit 
lime  scheme  for  Eq.  (B2)  would  yield  only  Z.  -I-  Z,j  „,(if.')  and  require  an  inversion  of  the 
operator  to  obtain  0'.  A belter  choice  is  an  implicit  time  scheme  which  utilizes  the  same  inver- 
sion. An  ideal  candidate  is  the  implicit  trapezoidal  scheme  described  by  Kurihara  (1965). 
Another  advantage  of  this  scheme  is  that  the  operator  inversion  required  for  Z.  | -f-  Z,  can  be 
used  to  compute  steady-stale  solutions  to  Eq.  (Bl)  without  wave  forcing  as  initial  conditions  for 
the  time-dependent  problems.  The  coupling  terms  between  the  wave  motion  and  mean  circula- 
tion can  be  made  second  order  accurate  in  time  by  computing  the  wave  field  at  times 
different  from  the  mean  flow. 


If  Ci-  L,(0)  and«j,„-  Lj,„(0) 


Then  Eqs.  (Bl)  and  (B2)  can  be  written,  respectively,  as: 


Ar 
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(B3) 


and 


4 » 4.  / 0 I 

SI  + S2.zn  + ^ 


‘frV  + -I-  T„-0“  + ^,0“  + 

or 


- y 


oz 


(B4) 


28 


r 1 

NRI  Ml  MOR  ANDl'M  RI  POHl  4013 
where  superscripts  denote  the  following  points  in  time: 

()*-() 

()  -()  A, 

1 

n--  ( „ 

Evaluation  of  the  spatial  operators  is  done  on  a staggered  mesh  spaced  at  5°  from  pole  to 
pole  with  the  closest  pole  points  at  87.5°,  for  ij/  and  85°,  for  li/'.  The  vertical  grid  spacing  is  0.37 
scale  heights.  The  generalized  sweepout  method  developed  by  Madala  (1978)  is  used  to  invert 
the  spatial,  elliptic  operators.  The  nonlinear  forcing  of  Eq.  (B3)  restricts  the  time  step  to  the 
equivalent  leapfrog  time  scheme  limit.  To  obtain  the  correct  group  velocity  for  planetary  waves 
with  the  implicit  trapezoidal  scheme  the  C.F.L.  condition  must  be  satisfied.  On  the  basis  of 
.several  numerical  experiments,  the  maximum  allowable  time  step  was  found  to  be  about  0.12 
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